A Cellular Automaton Model for Bi-Directional Traffic 
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We investigate a cellular automaton (CA) model of traffic 
on a bi-directional two-lane road. Our model is an extension 
of the one-lane CA model of jl| , modified to account for inter- 
actions mediated by passing, and for a distribution of vehicle 
speeds. We chose values for the various parameters to approx- 
imate the behavior of real traffic. The density-flow diagram 
for the bi-directional model is compared to that of a one-lane 
model, showing the interaction of the two lanes. Results were 
also compared to experimental data, showing close agreement. 
This model helps bridge the gap between simplified cellular 
automata models and the complexity of real-world traffic. 



I. INTRODUCTION 

A variety of cellular automaton models give a discrete 
approximation to traffic flow. Nagel and Schreckenberg's 
|jj one-lane model accurately describes some of the main 
features of real traffic. This lead to multi-lane S and 
two-dimensional (e.g. HQ) CA models, which extended 
the range of traffic phenomena which can be treated by 
cellular automata. This paper continues in this line by 
introducing a CA model in which vehicles move on two 
lanes, in opposite directions. We find that even com- 
plicated interactions, such as occur during passing into 
oncoming traffic, can be represented in a CA model. This 
supports the view that cellular automata provide a real 
alternative to differential equations for the modeling of 
traffic flow. 

Multi-lane models provide rules for lane changes, and 
we observe non-trivial interactions between the lanes. 
With two lanes in opposite directions, these interactions 
depend strongly on the relationship between the densities 
on the two lanes. When either or both of the densities 
are large, interactions can be ignored since passing is im- 
possible; our model is equivalent to two copies of the 
standard one-lane model. If there are no vehicles in the 
on-coming lane, then our model behaves like an asym- 
metric model with two lanes in the same direction. If 
there are only a few vehicles in a given lane and an in- 
termediate number in the oncoming lane, then vehicles 
in the oncoming lane will pass, and slow progress on the 
given lane. 

The one-lane CA model of PJ assumes that all vehicles 
have the same maximum velocity. Indeed, even a multi- 
velocity model quickly settles down to a state where all 
vehicles move with the lowest speed; the model unre- 
alistically creates platoons of vehicles, each following a 



slow car. As will be seen below, a distribution of speeds 
is especially important for realism in the bi-directional 
model. 



II. EXISTING MODELS 

a. One Lane Model In the one-lane CA model of Q| 
each site may be empty or occupied by a car with an in- 
teger velocity v 6 {0 . . . v max }. v max = 5 or greater gives 
good agreement with physical experiments. The variable 
gap gives the number of unoccupied sites in front of a 
vehicle. Pdecei is the probability to randomly decelerate, 
and rand is a random number between and 1. One 
iteration consists of the three following sequential steps, 
which are applied in parallel to all cars: 

1. Acceleration of free vehicles: IF (v < v max ) THEN 
v — v + 1 

2. Slowing down due to other cars: IF (v > gap) 
THEN v = gap 

3. stochastic driver behavior: IF (v > 0) AND ( rand 
< Pdecel ) THEN v = v - 1 

These simple conditions already give realistic results. 

b. Two-Lane Uni- directional Model The two-lane 
uni-directional model is built from two parallel single- 
lane models Four additional conditions govern the 
exchange of vehicles between lanes. First the vehi- 
cles change lanes, then the one-lane algorithm is ap- 
plied. This model introduces several new variables: 
gap S ame{i), gap pp(i)' the number of unoccupied sites in 
front of vehicle i on the (same, opposite) lane, respec- 
tively. gapbehind(i) the number of unoccupied sites behind 
the vehicle, on the opposite lane; l sa me,lopp, hack- the 
minimum free distance needed for a pass, ahead on same 
lane, ahead and behind on the opposite lane; p c hange'- 
the probability to change lanes. The added rule is as 
following: 



4. IF (gap same {i) < heme) AND (gap opp (i) > l opp ) 
AND (gap behmd (i) > hack) AND (rand < p c hange) 
THEN change lane 

This model has a symmetric and an asymmetric ver- 
sion. In the asymmetric version, there is no passing on 
the right. 
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III. TWO-LANE BI-DIRECTIONAL MODEL 

The bi-directional model has several types. Either- 
passing is allowed on both lanes, only one or the other 
of the lanes, or no passing is allowed. Mixtures of these 
types are possible; one type may be applied on a given 
section of roadway, and another type on an adjacent sec- 
tion. 

Consider a type in which passing is allowed. An algo- 
rithm which controls passing must account for a number 
of circumstances. First, a vehicle must not decelerate 
while passing; p — during the pass. Moreover, if an 
on-coming car is seen, the passing car must immediately 
return to its own lane. In reality vehicles generally do 
not attempt to pass unless the pass can be completed; a 
model should reflect this. We let each vehicle measure 
the " local density" : the density of cars in front of the 
vehicle it would like to pass. If the local density is suffi- 
ciently low, the vehicle has a good chance of completing 
a pass, and we allow it to try. Few passes will occur 
when the global density is high, even if the density is 
low on the on-coming lane. The lanes become effectively 
de-coupled. 

Three kinds of traffic jams can appear on a bi- 
directional road. The most common is a start-stop wave 
on one of the two lanes. Rarely, a jam is caused by an 
audacious driver tries to pass though there is no space 
to return to the home lane. Finally, a jam may occur 
when each of an adjacent pair of drivers, one on each 
lane, tries to pass simultaneously. Unless the symmetry 
between lanes is broken (see below), this "super jam" 
halts all traffic on both lanes. 

Our CA operates on a lattice with two lanes in op- 
posite directions (see figure 1). Each site has a state 
v G { — (v max + 1) . . . v max + 1}, where represents the 
absence of a vehicle, ±1 a stopped vehicle, ±2 a vehi- 
cle moving with speed 1 in the positive, resp. negative 
direction and so on. 

Vehicle movement is calculated in a two-step process, 
following First vehicles change lanes, then they ad- 
vance. 

In addition to the functions and variables introduced 
above, we will need the following: v same ,v opp : the ve- 
locity of the vehicle ahead on the (same, opposite) lane. 
H: true iff the vehicle is on its home lane; oncoming: 
true iff sign(v sam e) ^ sign(v); / pass : if gap same < Ipass 
AND H then a pass may be attempted. Iback- the dis- 
tance a driver looks back for obstacles on the passing 

lane. I security- it gap same < Security AND not{H) then 

the vehicle returns immediately to its home lane. D^: lo- 
cal density: the fraction of the Idensity — 2 x v max + 1 sites 
in front of the given vehicle which are occupied; Du m i t : 
the maximum local density for a safe pass. 

Spacel: true if (gap S ame < Ipass) AND (gap opp > 

^security) AND (gapbehind ^ ^back) ■ 

Space2: true if (gap opp > I security) AND (gap be hind > 

Iback) ■ 



In our simulations we use the following values through- 
it. Ipass ^3 Iback V nlaX7 ^security 2 X f maa; ~\~ 1, Dl 
T, — ! Pchange = 0.7, and Pdecel =0.5. 
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FIG. 1. bi-directional model 

Lane changes are determined by the following steps: 

1. IF (H AND spacel AND (D L < D Umit ) AND 
(rand < p c hange)) THEN change lane 

2. IF (not(H) AND {(gap same < hecunty) OR 
Space2)) THEN change lane 

The first condition affects vehicles on their home lane. 
If a vehicle is in front of them, of the same sign, and at 
a distance less than l paS s, then they would like to pass 
that vehicle. However, a pass will only be initiated if 
there is room far enough ahead on the passing lane, and 
the number of cars in front of the vehicle it would like to 
pass is small. Passing occurs randomly, even if all these 
conditions are met, the probability of changing lanes is 
denoted p c hange ■ The second condition concerns vehicles 
in the midst of passing. They return to their home lane 
if forced to by an oncoming vehicle, or if there is space 
enough on the home lane that they can return without 
braking. 

Forward motion of a vehicle is determined as follows: 

1. IF (M ^ v max ) THEN v = v + sign{v) 

2. IF ((oncoming) AND (gap sa me < (2 x v max - 1))) 
THEN v = [gaps am e/2\ 

3. IF ( (not(oncoming)) AND (|u| > gap sam e)) THEN 
v = sign(v) x gap sa me 

4. IF (( H AND (\v\ > 1) AND (rand < p deC ei) AND 
not(oncoming)) THEN v = V — sign(v) 

5. IF ( H AND (oncoming) AND (\v\ > 1) ) v = 
v — sign(v) 

These rules (1) accelerate the vehicle to maximum ve- 
locity, (2) rapidly decelerate the vehicle if there is an 
oncoming car too close, (3) decelerate the vehicle if it 
is closing in on another, both in their home lane, and 
(4) randomly decelerates the vehicle if it is on its home 
lane; if it is passing, it never decelerates randomly. Fi- 
nally, (5) breaks the symmetry between the lanes, and 
thus prevents the emergence of a super jam. 
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IV. RESULTS 

The usual way to represent the behavior of a CA is with 
a space-time diagram, where space is represented on the 
horizontal axis, and time on the vertical axis, with time 
proceeding downward. A typical space-time diagram for 
this model is shown in figure 0. 

Space 



This experiment was performed with a value of Pdecei of 



passing zone 



no passing zone 




FIG. 2. passing and no-passing zones. 
A cell-wide section of roadway containing only one vehicle 
moving (rightward resp. leftward) on its home lane is labeled 
red resp. green). Other colors give other possible situations, 
such as a single vehicle passing, going rightward (violet). On 
the left-going lane the density is near critical for jams (0.1) 
and the density is low (0.01) on the right-going lane (500 cells, 
periodic boundary conditions). The route is divided into a 
passing zone (left half of the figure) and a no-passing zone 
(right half of the figure). Jams occur much more frequently 
in the no-passing zone. All vehicles have the same maximum 
velocity of 5. 

This figure illustrates that passing can dramatically 
fluidity traffic. The start-stop waves seen in the right 
half of the figure disappear in the left half. Here the 
densities have been chosen to maximize the effects caused 
by interactions between two lanes in opposite directions. 

We have explored the entire range of densities on the 
two lanes, as shown in figure IV. Here the difference 



between the flow on the home lane in two-lane model 
with passing on both lanes is compared to the flow in 
a one-lane model (A/). When the density on either or 
both lanes is large then there is little difference between 
the two-lane and one-lane models. When the density on 
the passing lane is small ( < 0.1), then the flow on the 
home lane can be much greater than in a one-lane model, 
blue). Maximum improvement occurs near density on 
the passing lane. If the density on the home lane is small 
( < 0.25) then the flow may be lower than in the corre- 
sponding one-lane model since when oncoming cars pass 
other oncoming cars they can impede traffic on the home 
lane. In an asymmetric model in which only vehicles on 
the home lane can pass, this slowing effect disappears. 



0.5 and Pchange of 0.5 




FIG. 3. 3D plot of flow on the home lane. 
The height of the surface gives the difference between flow 
on the home lane in the bi-directional model and flow in a 
one-lane model as a function of the densities on the home 
and passing lanes of the bi-directional model. 

Our bi-directional model allows us to treat a host of 
new phenomena. For instance, the multi-speed variant 
(figure |^) can produce a complicated knot of interac- 
tions resulting from a slow car leading a group of oth- 
ers through a no-passing zone. While in the no-passing 
zone the faster cars bunch up behind the slow car, form- 
ing a platoon. When the platoon reaches the end of the 
no-passing zone the platoon disperses. 

Space 



passing zone no passing zone passing zone 



FIG. 4. A space-time diagram generated by a multi-speed 
model 

Each car has a maximum velocity, uniformly dis- 
tributed between 2 and 5. In the center half of the figure 
no-passing is allowed, while on the two sides passing is 
allowed. Slower vehicles trace steeper lines. Note in par- 
ticular the large knot of cars near the center of the figure 
formed by a slow (left-going) car entering a no-passing 
zone. As the knot passes the leftmost boundary of the 
no-passing zone, it disperses. 
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V. COMPARISON WITH EXPERIMENT 

In figure [^, we compare simulation results (left panel) 
with the measurements of M. Van Aerde and S.Yagar 
H (right panel) . These experimental data are supported 
by traffic volume data from Australia and Canada ||. 
Van Aerde and Yagar measured the maximum flow which 
can be obtained on the home lane, as a function of the 
on-coming flow on the passing lane. We made compara- 
ble measurements on our simulations, with a variety of 
parameter setting, so as to approximate the parameter 
values implicit in the experimental data. 
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accounted for by a few simple rules. By small extension 
of previous work, we have greatly increased the realism 
which can be expected from cellular automaton models 
of vehicular traffic. Traffic in the two-lane bi-directional 
model approximates real traffic on both a microscopic 
and macroscopic level. Indeed, our results suggest revi- 
sion of a standard text in highway studies. 

Though systematic exploration of the parameter space 
should surely be undertaken, our model is robust to small 
changes in parameter value. Throughout the experiments 
described in this paper, we have fixed the values of most 
parameters at reasonable values, and focused on the two 
most important variables: the density of vehicles on each 
lane. We have localized the regions where lane interac- 
tions are strongest, and explored how these interactions 
occur. Due to the high relative speeds involved, these 
interactions have practical importance in traffic safety. 
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FIG. 5. Comparison of simulation results with experimen- 
tal data. 

left panel: The results of simulations under three conditions, 

a)pd ece i = 0.5, Pchange = 0.5, b)pdecsl = 0.25, Pchange = 0, 

and c)pdecei = 0.25, Pchange = 0.3. right panel: corresponding 
physical measurements, modified from figure 5 of [5]. 

In the simulations, we obtain maximum total flow 
when vehicles are distributed uniformly, that is, half in 
each direction. The simulation curves are obtained by 
varying the total number of vehicles. 

In both physical experiment and computer simulation, 
when densities are high, flows are not affected by pass- 
ing, resulting in a corner in the graph of maximum flow. 
When one of the lane densities is fixed to a small enough 
value, the maximum flow increases on the other lane due 
to passing. Note that while our simulations are in close 
qualitative agreement with these experimental results of 
H , these same results are not in agreement with the re- 
sults reported in the previous jf| and current editions 
of the influential Highway Capacity Manual (HCM) . The 
previous HCM || presents a straight line of slope -1, 
while the current HCM presents a smoother curve, 
which suggests that flows are correlated at high densities. 
Our simulation, built from close modeling of microscopic 
interactions, gives theoretical support for the view that 
the HCM's treatment of bi-directional traffic should be 
revised in favor of the results of Van Aerde and Yagar. 



[1] K. Nagel and M. Schreckenberg, J. Phys. I France 2, 2221 
(1992). 

[2] M. Rickert, K. Nagel, M. Schreckenberg, and A. Latour, 

Physica A 231, 534 (1996). 
[3] J. A. Cuesta, F. C. Martinez, J. M. Molera, and A. 

Sanchez, Phys. Rev. E 48(6), R4175 (1993). 
[4] T. Nagatani, Phys. Rev. E 48, 3290 (1993). 
[5] S. Yagar, Australian Road Research 13(1), 3 (1983). 
[6] Highway Capacity Manual 1965. HRB Special report 87. 

HRB, National Research Coucil, Washington, D.C. 397 pp 

(1965). 

[7] Highway Capacity Manual 1994, Special Report 209, 
Third Edition, TRB, National Research Coucil, Washing- 
ton, D.C. chapter 8 p6 (1994). 



VI. DISCUSSION 



We have shown that complicated interactions between 
vehicles traveling on two-lane rural roads can be largely 
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